Semiclassical Theory of Time-Reversal Focusing 
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Time reversal mirrors have been successfully implemented for various kinds of waves propagating 
in complex media. In particular, acoustic waves in chaotic cavities exhibit a refocalization that is 
extremely robust against external perturbations or the partial use of the available information. We 
develop a semiclassical approach in order to quantitatively describe the refocusing signal resulting 
from an initially localized wave-packet. The time-dependent reconstructed signal grows linearly 
with the temporal window of injection, in agreement with the acoustic experiments, and reaches the 
same spatial extension of the original wave-packet. We explain the crucial role played by the chaotic 
dynamics for the reconstruction of the signal and its stability against external perturbations. 
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The concept of time reversal has captured the imagi- 
nation of physicists for more than a century, leading to 
a vast theoretical oeuvre, sempiternal discussions, and 
a few concrete experimental realizations. Among them, 
the works on spin echoes have been of paramount impor- 
tance concerning the limits in the reconstruction of an 
initially prepared quantum state [l|, |3] ■ The time rever- 
sal of acoustic waves in a non-homogeneous medium was 
another experimental deed showing that an initially lo- 
calized pulse can be accurately reconstructed by an array 
of receiver-emitter transducers that re-inject the recorded 
signal ^]. Re-focusing of elastic, as well as electromag- 
netic, waves has been later achieved [J, S Ifil- These ex- 
periments provoke a natural surprise while yielding re- 
constructions, that albeit not perfect, are highly faithful. 
The relevant questions that arise when trying to under- 
stand these physical realizations of time reversal are re- 
lated with how good a reconstruction can be achieved and 
which are the limits set by interactions with the environ- 
ment and unavoidable errors in the reversal protocol. 

In the spin echo experiments the complexity of the 
physical system has emerged as a critical component, 
and the term of Loschmidt echo (LE) has been coined 
to describe setups where many-body physics or chaotic 
dynamics are relevant Q . In particular, the decay of the 
LE with the reversal time has been shown to depend on 
the underlying classical dynamics [3, |§, [l0|: classically 
chaotic systems exhibit a decay rate which, for large per- 
turbations, is bounded by the Lyapunov exponent char- 
acterizing the dynamics. 

In the time-reversal mirror (TRM) procedure the play 
back signal builds up in the region of the original ex- 
citation, in the form of a reversed wave amplitude [3|. 
Thus, the TRM can be viewed as the wave version of 
the LE. A salient feature of the TRM experiments is 
that, even though reversal is not perfect ll|, the refo- 
cusing improves when the wave-propagation occurs in a 
disordered medium or in a chaotic cavity, as compared 
with the homogeneous or integrable case. Remarkably, 



a single transducer is enough in the case of a chaotic 
cavity 4j . The asymptotic analysis of the Wigner trans- 
form of wave fields in the high-frequency limit has been 
used to understand how multiple scattering enhances the 
spatial resolution of the refocused signal [12]. Diagram- 
matic perturbation theory has been able to account for 
the symmetry-induced interference enhancements in the 
refocalization observed in disordered media [13| . The re- 
focusing experiments in chaotic cavities have been con- 
fronted with numerical simulations [J], as well as ergod- 
icity and control theory [14l|. The contrasting stability 
properties of TRM with wave and particle propagation 
through a multiple scattering medium has been discussed 
in Ref. [l5|. 

In this work we develop a semiclassical approach for 
TRM in chaotic cavities and quantify the quality of the 
reconstructed signal in terms of temporal and spatial dis- 
persions, as well as possible environmental influences. We 
demonstrate the crucial role played by the underlying 
classical dynamics and validate our analytical results by 
confronting them to numerical simulations. 

A high-frequency signal emitted at i = at a posi- 
tion ro inside the cavity can be interpreted within the 
ray picture as an initial wave-packet centered at Tq that 
evolves and is recorded by a receiver (or an array of re- 
ceivers) at position(s) r^ {i = 1,2,..., TV) for times in 
the interval (ii, ^2)- After a waiting time tw > ^2 the re- 
emission of the time-reversed signal is performed in the 
interval {t'2 — 2tw — h^t'i = 2tw — ti)- The refocusing 
is expected at 2t\\r, that is redefined as the time origin 
for refocusing [J| (see Fig. [TJa)). The signal that can be 



detected in a point r at a time t is |16| 



•^Po(l",t) = 



A^ 



dTG(r,r„i + T) 



X /dr'G*(r,;,r',T)V';„(r') 



(1) 



The propagator G{v,ri,t + r) corresponds to the re- 
emitted signal, which is obtained by time-reversing 



the evolution of the initial state with the propagator 
G(ri,r',r). We do not write the initial temporal argu- 
ments of the propagators, as they are taken to be 0. We 
work in two dimensions and choose as an initial state a 
Gaussian wave-packet 

1 



V'po(r') 



/TTcr 



■ exp 



iPo • (r' - ro) 






(2) 



centered around ro and with dispersion o. The initial 
momentum po sets the energy and direction of the origi- 
nal signal. The choice of a quantum formalism to repre- 
sent the ray picture is motivated by convenience, as we 
are leaving aside the delicate issue concerning a quantal 
recording-emission process. The quantum formalism is 
suited to work with the semiclassical propagator [l7| 



G(r",r',r)== 
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(3) 



s(r', 



r) 



given as a sum over classical trajectories s traveling in a 
time T between the two extreme points. We note Ss = 
Sg (r" , r' , r) the action integral along the path, /i^ the 
Maslov index, and C, = \det{-dSs/dr'dr")\. 

Performing the r'-integral of Eq. ([T]) by stationary- 
phase (see Ref. J^j for the precise conditions of valid- 
ity of such an approximation in a chaotic cavity) we can 
write, for a single transducer, 



J^po(r,t) 
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dr 
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where p^ is the initial momentum of the trajectory s. 

We are interested in times t close to the refocusing 
one, and positions r near ro, therefore the dominating 
contribution comes from the diagonal approximation s' ~ 
s leading to a signal given by 



•^Po(i",t) 



'27r3/2?i2 



dr Fpo(ro,rj,r,r), (5a) 



i^p„(r',r",r,r)= }_^ C, /p„(r', r",r,p'), (5b) 

s(r' ,r" ,r) 



/p„(r',r",r,p')=exp[ip'.(r-r') 



-j-Est 



i-(p'-Po) 
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(5c) 



Es is the energy at which the trajectory s is traveled. 
In billiards the magnitude of the momentum modifies 
the traveling time but does not affect the path. Not- 
ing s(ro,ri) the geometrical support of the trajectory 
s(ro,ri,r) with length Lg, we have Eg — p^/2m and 
Ps = (rriLg/T), where m is the mass of the particle. 

In order to present the calculation in its simplest form 
we start by setting po = and the optimal condi- 
tions t = and r = ro, which yield (from the last 



term in the exponent of ([5c)) ') the maximum refocusing 
•?^max = •?^o(ro, 0). In a fully chaotic system, Cg scales as 
exp[— Ast], where As is the largest Lyapunov exponent. 
Assuming further a uniformly hyperbolic dynamics [l8| 
and using that in a billiard XgT = XLg (with A as an 
inverse length), we write Cg = 2m'^XLs/T^ ex.p[—XLs]. 
Noting V^ = 2Tt? / rri^ cr"^ we have 
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The sum over the transient orbits s can be converted 
into an integral over trajectory lengths by introducing the 
density dN{L)/dL = 7r/(Ayt) exp (A£) (A stands for the 



area of the chaotic cavity) 



Denoting Lj = Vtj 



{j = 1,2) and Ld the length of the shortest trajectory 
linking ro and r; we have 



-^max — 7 { (^2 ^ ^l) 

aA 



id/i2 



dl [1 - erf {I)] 



+ii / dl[l~ erf (/)] ) , (7) 

"'id/L2 J 

where erf (x) stands for the error function. The assump- 
tions made on Cg and dN{L)/dL are valid for lengths a 
few times larger than L^. However, our approximation is 
appropriate since we assume that we start recording at 
times ti large enough for the typical contributing trajec- 
tories to feel the chaotic nature of the dynamics. That 
is, we work under the hypothesis L^ <^ Li < L2, that 
also allows to neglect the last integral, leading to 



J'n 



/naA 



AT . 



(8) 



The scaling of the refocused signal with the injection in- 
terval AT = ^2 — ^1 is a quite natural result, experimen- 
tally observed in Ref. [4i], while the scaling with A has 
not been systematically tested so far. In the case where 
there is an array with N transducers, we simply have to 
multiply the above result by iV, but the surprising fact 
that just one detector is enough stems from Eq. ([8]). In 
order to further quantify faithfulness of the time-reversal 
process we need to evaluate the temporal and spatial ex- 
tents of the reconstructed signal. 

For times t close to the refocusing one we have to 
consider the phase mL'jt/{2T'^h) in Eq. (|5c|) . Therefore, 
To (ro , t) follows the same expression than J^max in Eq. ([6]) 
by the change of l/V^ into l/V^ + imt/2h. The error 
functions resulting from the r-integral have now to be 
extended into the complex plane, and after a straightfor- 
ward algebra we find 



To{ro,t) = 



T„ 



l + iVV(\/2cr) 



(9) 



The reader can imagine how the general calculation 
goes when we treat simultaneously i ^ 0, r ^ Tq and 
Po ^ 0. Instead of presenting such calculation, which 
results in the faithful reconstruction of the initial wave- 
packet, we look at the problem from a different perspec- 
tive and introduce the ergodicity hypothesis in order to 
treat the general case. The ergodic approach not only 
provides a second, and more economical, way of obtaining 
the general result without using detailed knowledge of the 
dynamics, but also sheds some light into the necessary 
conditions for achieving the refocalization condition. The 
basics of the ergodic approach is to calculate quantities 
like Fjvjj(r',r",r,T) of Eq. ((5b| by averages over phase- 
space [21|. Calling Vr — rT-(r',p') and Pr = PT(r',p') 
the position and momentum, respectively, at time r of a 
particle starting at (r', p') at time 0, we have 

i^p„ (r', r", r, r) - 1 dp'dp"^ (r. - r") <5 (p. - p") 

x/p„(r',r",r,p'). (10) 

The double delta- function represents the distribution 
of classical trajectories. An average over small ranges of 
initial and final conditions gives a smooth distribution 
which describes the evolution in a statistical sense. For 
sufficiently long times such a distribution is r indepen- 
dent, and uniformly distributed on the hyper-surface of 
constant energy (which for two dimensional billiards has 
a volume Q. — 2TTmA in phase-space) . We therefore have 



Excitation Recording 



Fp„(r',r",r,r)=|dp'dp"<5(H:^) 

= iy'dp'/p„(r',r",r,p'). 



/p„(r',r",r,p') 
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(11) 



Applying this general procedure to the function /p^ of 
Eq. ([Sc]) we have 



•^Po(i-, t) =^^J^2j J dp' exp [ip' • (r - ro) 



^2m~ 



-t 



2h- 



■(p'-Po) 



(12) 



since the integral over r is now trivial. Performing the 
Gaussian integral over p' we obtain a wave-packet that 
refocalizes with the same shape of the original one, but 
with momentum — po. The magnitude of the signal close 
to the maximum refocalization condition is given by 



l-^Po(r,t)l = 



T„ 



^1 + {Vt/V2ay 



: exp 



{^-^0 + ^t) 

2^2 + (Vt)2 



(13) 

Numerical calculations of time-reversal imaging have 
been performed in Ref. [J| for a two-dimensional cavity 
with the shape of a sliced disk. The signal reconstruc- 
tion could be visualized and a qualitative agreement with 
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FIG. 1: (color online) . (a) TRM sequence, (b) Reconstructed 
signal scaled with A/ AT at the emission point ro, close to the 
refocusing time f = 0, for the shown billiard. The thick solid is 
the semiclassical prediction (Eq. (|13p ). Numerical simulations 
for various AT and A: 5000 and 150x300 (blue dotted), 10000 
and 150 x 300 (green dashed), and 5000 and 300 x 600 (red 
dash-dotted). Right inset: reconstructed signal at f = close 
to ro from the semiclassical prediction (solid) and simulation 
(red dash-dotted) with AT = 5000 and ^ = 300 x 600. 



the experimental results was found. Since we dispose 
now of a quantitative semiclassical theory of refocusing 
it is important to test our predictions in a stadium bil- 
liard, which is a paradigm of classical chaotic dynamics. 
We calculate the evolution of the wave-packet through 
a second order Trotter-Suzuki algorithm for a discrete 
Schrodinger equation. Lattice effects are minimized by 
considering ao ^ Ab <C cr ^ Lb, where ag is the lattice 
constant, Ab the de Broglie wave-length associated with 
Po, and ifc the size of the billiard. We assume that at 
injection time all the original signal has already decayed. 
In Fig. [IJb) we show the numerical results for the time 
dependence on the reconstructed signal at rg. The nor- 
malized J- A/ AT are well described by the semiclassical 
prediction (thick solid) confirming the scaling with AT 
and A of Eq. ^ . The normalizing factor for the numer- 
ical results is approximately 1.4 times the semiclassical 
one. Such a difference may be due to our discretiza- 
tion of the quantum problem as well as the difficulties of 
the diagonal approximation to recover exact numerical 
values. The signal-to-noise ratio does not change appre- 
ciably when the recording time is doubled, while it is 
improved by increasing A. In the right inset we show 
the spatial reconstruction of the wave-packet around ro 
for the refocusing time t = 0. We see that the semiclas- 
sical prediction (thick solid) provides the proper scaling 
behavior and, up to the normalization factor, a quanti- 
tative description of the TRM results. 



The numerical implementation of TRM for integrable 
geometries results in a refocusing that strongly depends 
on the position of the transducers and with a signal 
hardly distinguishable from the background. The semi- 
classical approach allows to understand this important 
difference between chaotic and integrable systems. In 
the former the exponential proliferation of trajectories 
allows to encode the information at all times, while in 
the latter the registered signal will be strongly depen- 
dent on whether or not the source and the transducer 
are connected by a stable trajectory. 

Experimentally, the TRM procedure has been shown 
to be robust against local and global perturbations in- 
troduced between the recording and injection phases [3|. 
Even in the absence of these perturbations, it is natu- 
ral to expect that in any TRM setup the environment 
acting during the re-emission is slightly modified respect 
to that of the recording phase. In the same spirit of 
the LE studies, we can model this situation by assuming 
that in the recording process we have a Hamiltonian H 
that determines G(ri,r',r) in Eq. ([T]), while a modified 
G(r, Yi,t + T) is governed by the slightly different Hamil- 
tonian H . For the LE the details of the perturbation 
H — H are not important, and its effect can be accounted, 
after some averaging, by affecting the contribution of 
each trajectory s by an additional factor exp[— Ls/2/], 
where I is an effective mean-free-path associated to the 
perturbation. In general I depends on the velocity of the 
particle, i.e. for perturbations modeled by an auxiliary 
impurity potential characterized by a strength 7 we have 
l/l = ^t'^ /Li f!!]. Including this r-dependent phase pre- 
vents us of using the ergodic approach, but working along 
the lines of Eqs. (IT])-® we obtain a maximum refocusing 
for a non-static environment given by 
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exp 
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exp 



1L2 



(14) 



which reduces to JF^^x (1 — {h +ii)/4T) for t2 ^ f; 
to .F,„axexp[-(i2 +ii)/4f] for ^2-^1 < f; and to 
CT/{y/TTaA) exp [—c'ti/2f] for ^2 >'r. The numerical con- 
stants c = 2.94 and c' = 0.46 are given by rational- 
argument values of the F function, and the characteristic 
time is defined as f = V/(2y^7). The proportionality 
of the reconstructed signal with the injection interval is 
clearly lost in the limit of ^2 > "? since the perturbation 
renders ineffective the longest recording times. From the 
relevant limits that we have singled out, the second one 
is the most important for current experiments. It has a 
Fermi-Golden-Rule structure that can be obtained under 
very general considerations. For perturbations where the 
effective elastic mean-free-path / increases with the Lya- 
punov exponent of the unperturbed system (i.e. mass 



distortion in a Lorentz gas 9|) the characteristic time f 
increases with the chaoticity of the system. Similarly to 
the Fermi- Golden-Rule regime of the LE, such a motional 
narrowing effect translates into larger stability, and im- 
proved refocalization, for the more chaotic systems, in 
agreement with the experimental findings [3[ . 

In summary, we have described the refocalization sig- 
nal for the time reversal mirror procedure through the 
semiclassical approximation. The chaotic nature of the 
underlying classical dynamics appears as a key ingredi- 
ent to ensure the stability of the refocalization towards 
perturbations and the proportionality of the recovered 
signal with the injection time. 
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